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ABSTRACT 


Context. A long-held vision has been to realize diffraction-limited optical aperture synthesis over kilometer baselines. This will en¬ 
able imaging of stellar surfaces and their environments, and reveal interacting gas flows in binary systems. An opportunity is now 
opening up with the large telescope arrays primarily erected for measuring Cherenkov light in air induced by gamma rays. With 
suitable software, such telescopes could be electronically connected and also used for intensity interferometry. Second-order spatial 
coherence of light is obtained by cross correlating intensity fluctuations measured in different pairs of telescopes. With no optical links 
between them, the error budget is set by the electronic time resolution of a few nanoseconds. Corresponding light-travel distances are 
approximately one meter, making the method practically immune to atmospheric turbulence or optical imperfections, permitting both 
very long baselines and observing at short optical wavelengths. 

Aims. Previous theoretical modeling has shown that full images should be possible to retrieve from observations with such telescope 
arrays. This project aims at verifying diffraction-limited imaging experimentally with groups of detached and independent optical 
telescopes. 

Methods. In a large optics laboratory, artificial stars (single and double, round and elliptic) were observed by an array of small tele¬ 
scopes. Using high-speed photon-counting solid-state detectors and real-time electronics, intensity fluctuations were cross-correlated 
over up to 180 baselines between pairs of telescopes, producing coherence maps across the interferometric Fourier-transform plane. 
Results. These interferometric measurements were used to extract parameters about the simulated stars, and to reconstruct their two- 
dimensional images. As far as we are aware, these are the first diffraction-limited images obtained from an optical array only linked 
by electronic software, with no optical connections between the telescopes. 

Conclusions. These experiments serve to verify the concepts for long-baseline aperture synthesis in the optical (somewhat analogous 
to radio interferometry) and to optimize the instrumentation and observing procedures for future observations with large arrays of 
Cherenkov telescopes, aiming at order-of-magnitude improvements of the angular resolution in optical astronomy. 

Key words. Intensity interferometry - Hanbury Brown - Twiss - Aperture synthesis - Second-order coherence - Stellar surface 
imaging - Exoplanet imaging - Cherenkov telescopes 


1. Introduction 

Highest angular resolution currently realized in optical astron¬ 
omy is offered by amplitude (phase-) interferometers that com¬ 
bine light from telescopes separated by baselines of up to a few 
hundred meters. Tantalizing results show how stellar disks start 
to become resolved, revealing stars as a diversity of individual 
objects, although so far feasible only for a few of the largest 
ones. Several concepts have been proposed to extend such facili¬ 
ties to scales of a kilometer or more, but their realization remains 
challenging either on the ground or in space. Limiting parame¬ 
ters include the requirement of optical and atmospheric stability 
to a fraction of an optical wavelength, and also the need to cover 
many interferometric baselines, given that optical light cannot be 
copied with retained phase, but has to be split up (and diluted) by 
beamsplitters to achieve interference among multiple telescope 


pair s. Overviews of thos e techniques are given by Millour| ( 2008 ) 
and |Quirrenbach| { 2009|) and in more detail by Monnier ( 2003 ) 
a nd|Monnier & All^ ( |2013|), while re cent results are highlighted 
in |Berger et al. ( |20T2 ) and |van Belle ( 2012| ). 


Bright stars typically subtend diameters of only a few mil- 
liarcseconds (mas), requiring interferometry over many hun¬ 
dreds of meters or some kilometer to enable surface imaging. 
Telescope complexes of such an extent are now being con¬ 
structed on the ground as arrays of air Cherenkov telescopes. 
These are erected to measure brief flashes of visual Cherenkov 
light in air induced by energetic gamma rays but their optome¬ 
chanical performance is fully adequate for intensity interfer- 
ometry also, a technique already pioneered long ago ( |Hanbury| 
Brown[ |1974| ). Because it is essentially insensitive to atmo¬ 
spheric turbulence, this method permits both very long baselines 
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and observing at short optical wavelengths; only the lack of suit¬ 
ably large and well-distributed telescopes seem to have caused 
the method not to be pursued recently in astronomy. 

Currently, the largest such project is CTA, the Cherenkov 
Telescope Array, foreseen to have around 100 telescopes spread 
over a few square kilometers ( [Acharya et al.|Pol'S}|CTA|[2015| ). 
Among its envisioned uses is also that of an intensity interferom¬ 
eter, correlating intensity fluctuations between numerous tele¬ 
scope pairs across different baselines whil e applying timeshifts 
in sof tware to track sources across the sky ( [Dravins et al.[ |2012| 
2013| ). Since the signals are copied electronically (with no op¬ 


tical connection between telescopes), there is no loss of signal 
when forming additional baselines between any pairs of the nu¬ 
merous telescopes. 

Theoretical and numerical simulations of such observations 
have been made for several types of obj ects ([Nunez et al.j 


2012a|b[[Rou et al.|[2013[[Dolne et al.[[2014[ ). Here, we report on 


their experimental verification. A laboratory setup was prepared 
with artificial optical sources (stars), observed by an array of 
small telescopes equipped with nanosecond-resolving photon¬ 
counting solid-state detectors. The photon streams (reaching 
MHz levels, as for realistic telescope operations) were fed into 
digital firmware correlators in real time (with nanosecond res¬ 
olution, as required for actual telescope operations), comput¬ 
ing cross correlations between the ‘random’ intensity fluctua¬ 
tions simultaneously measured in many different pairs of tele¬ 
scopes. The degree of mutual correlation for any given base¬ 
line provides a measure of the second-order spatial coherence of 
the source at the corresponding spatial frequency (and thus the 
square of the ordinary first-order coherence and the square of the 
Fourier transform component of the sources brightness distribu¬ 
tion). Numerous telescope pairs of different baseline lengths and 
orientations fill in the interferometric (u,v)-plane, and the mea¬ 
sured data points build up a two-dimensional map of the second- 
order spatial coherence of the source, from which its image can 
be extracted. The experience from such laboratory experiments 
will serve to verify instrumental concepts and serve as input for 
specifying the procedures for future full-scale observations with 
large telescope arrays. 


2. Principles of intensity interferometry 

The technique of intensity interferometry was developed for 
the original purpose of measuring stellar size s, and a dedicated 
instrument was built at Narrabri, Australia ( [Hanbury Brown[ 
1974[ ). What is observed is the second-order coherence of light 
(i.e., that of intensity /, not of amplitude or phase) by mea¬ 
suring temporal correlations of arrival times between photons 
recorded in different telescopes. A star is observed with two 
separate telescopes, recording the (quantum-optically) random 
and very rapid intrinsic fluctuations in starlight (fully resolved 
on timescales equal to the temporal coherence of light, perhaps 
10“^^ s). With the telescopes close together, the fluctuations in 
both telescopes are correlated in time, but lose correlation for 
larger separations. 

The quantity measured is 

{hit) ■ hit)) = {hit)){hit))ii + \ri2\h (1) 

where () denotes temporal averaging and yu is the mutual co¬ 
herence function of light between locations 1 and 2, the quantity 
commonly measured in phase/amplitude interferometers (a re¬ 
lation that holds for each polarization). Compared to randomly 
fluctuating intensities, the correlation between intensities A and 


h is ‘enhanced’ by the coherence parameter and an intensity in¬ 
terferometer thus measures I 712 P with a certain electronic time 
resolution. This relation holds for ordinary thermal (chaotic, 
maximum-entropy, Gaussian) light, where the wave undergoes 
random phase jumps over timescales of its coherence time but 
not necessarily for light with different photon statistics (e.g., an 
ideal laser emits coherent light of constant intensity, without any 
phase jumps, and thus would not generate any reasonable signal 
in an intensity interferometer). 

The great observational advantage of intensity interferome¬ 
try is that it is practically insensitive to either atmospheric turbu¬ 
lence or to telescope optical imperfections, enabling very long 
baselines and observations at short optical wavelengths, even 
through large airmasses far away from the zenith. Since tele¬ 
scopes are only connected electronically, error budgets and re¬ 
quired precisions relate to electronic timescales of nanoseconds, 
and light-travel distances of tens of centimeters or meters rather 
than small fractions of an optical wavelength. For more detailed 
discussions of the principal workings of intensity interferometry, 
see mo nographs by [Labeyrie et al.[ ( [2006[ ), [Sah^ ( [2011[ ), or [Shih[ 
( 2011[ ), for example. 

There is a price to be paid, however, for this freedom from at¬ 
mospheric influences. Realistic time resolutions are much longer 
than typical optical coherence times for broad-band light, and 
any measured intensity-fluctuation signal is averaged over very 
many coherence times. Therefore, very good photon statistics 
are necessary for a reliable determination of the smeared-out sig¬ 
nal, requiring large photon fluxes (thus large telescopes); already 
the 6.5 m flux collectors in the original intensity interferometer 
at Narrabri were larger than any other optical telescope at that 
time. 


The technical specifications for air Cherenkov telescopes 
such as the forthcoming Cherenkov Telescope Array ( [Acharya[ 
et al.[[2013[[CTA[[2015[ ), are remarkably similar to the require¬ 


ments for intensity interferometry. In the original Narrabri in 
strument, the telescopes were moved on tracks during observa¬ 
tion to maintain their projected baseline; however electronic time 
delays can now be used instead to compensate for different ar¬ 
rival times of a wavefront to the different telescopes, removing 
the need for having them mechanically mobile. 

With perhaps 100 telescopes, the CTA will provide an un¬ 
precedented light-collecting area of some 10,000 m^ distributed 
over a few square kilometers. Of course, it will mainly be de¬ 
voted to its primary task of observing Cherenkov light in air; 
however several other applications have been envisioned, to be 
preferably carried out during nights with bright moonlight which 
- owing to the faintness of the Cherenkov light flashes - might 
preclude their efficient observation. Besides intensity interfer¬ 
ometry ( [Dravins et al.| 2013[ ), additional uses that have been 


suggested include searches for rapid astrophysical events (Deil 
[et al.[ [2009t [Hinton et al.| [2006t [Lacki| [2011[ ), observations of 
stellar occultations by distant Kuiper-belt objects ( [Lackij 2014 ), 
or as a terrestrial ground station for optical communication with 
distant spacecraft ( [Carrasco-Casado et al.[ [2013[ ). For practical 
operations, one proposed concept is to place auxiliary detec¬ 
tors on the outside of the Cherenkov camera cover lid, an ap- 
proach already rea lized on some telescopes ( [Deil et ak] 


200 ^ 


Hinton et al. 2006[ ), and one that should not interfere with regu¬ 


lar Cherenkov camera operations. 


This potential of using Cherenkov telescope arrays for opti¬ 
cal intensity interferometry has indeed been noticed by several 
( [LeBohec & Holder] [2006| [LeBohec et~ar|[2008[ ). If, for exam¬ 
ple, a baseline of 2 km could be utilized with CTA at A 350 nm. 
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resolutions would approach 30 yuas, which is an unparalleled spa - 
tial resolution in optical astronomy ( [Dravins et al.[ 2012| 2013). 


3. Verifying intensity interferometry 


Several studies have considered the observational challenges, 
such as the overall sensitivity and limiting stellar magnitudes, 
and those issues now appear to be consistently and generally un¬ 
derstood ([Dravins et al.||2012||2013t|LeBohec & Holderj [2006| 


[Trippe et al.[ 2014[ ). For example, current types of electronics 
with current Cherenkov telescope designs should, during one 
night of observing in one single wavelength band, permit mea¬ 
surements of hotter stars down to a visual magnitude of about 
mv=8, giving access to many thousands of sources. 

Comprehensive Monte Carlo simulations were done by [Rou[ 
[et al.[ ( MT3] ). They numerically simulated effects of various noise 
sources, of background light from the night sky, effects of large 
collecting areas (i.e., non-negligible in comparison to the coher¬ 
ence pattern), among others. For example, the effect of large-size 
telescopes is that the signal will correspond to a convolution of 
the degree of coherence \y\^ with each of the telescope light- 
collection area shapes. However, in real observing, there will 
also arise various experimental issues, not all of which are still 
understood well. Some were examined already long ago during 
the Narrabri interferometer operations, while others are specific 
to modem detectors and electronics (dead-time or saturation ef¬ 
fects in detectors and digital data handling). 


3.1. Experimental test observations 


Full-scale experiments in connecting major Cherenkov tele¬ 
scopes for intensity interferometry were made by [Dravins 


& LeBohec ( 2008[ ), using pairs of 12-m telescopes of the 


VERITAS array on Mt.Hopkins, Arizona, at baselines between 
34 and 109 m. In each telescope, starlight was recorded by one 
photomultiplier in the center of the regular Cherenkov cameras, 
its signal digitized, and sent via optical cables to the control 
building, where they were fed into a real-time digital correlator 
at continuous photon-count rates up to some 30 MHz. Although 
no astrophysical correlations were recorded, these experiments 
verified that there seemed to be no particular operational prob¬ 
lems in using large Cherenkov telescopes for intensity interfer¬ 
ometry. 

Significant efforts have been made at The University of Utah 
in developing experimental test facilities for intensity interfer¬ 
ometry, including 3-m Cherenkov telescopes at its ‘StarBase’ 
facility and laboratory setups ( [LeBohec et al.[[20TQ| ). For exam¬ 
ple, an analog correlator based upon a Field Programmable Gate 
Array (FPGA) was used for laboratory tests ( [Nunez[[2012[ ) with 
intensity fluctuations in starlight mimicked by a pseudo-thermal 
light source, which is produced by shining a laser through a ro¬ 
tating ground-glass plate, producing a time-variable speckle pat¬ 
tern ( Martienssen & Spiller[[1964] ). 

Also [Pellizzari et al.[ ( [2012[ ) used a rotating ground-glass 
plate to create speckle patterns in the observation plane with two 
laboratory telescopes. Cross correlating the light intensity be¬ 
tween these provides a measure of the speckle extent, hence the 
source size. These experiments have provided insights into de¬ 
tector properties and signal handling although this type of light 
sources do not have the small angular extents that would be ap¬ 
propriate for astronomical imaging. 


3.2. Photon-counting detectors 

Air Cherenkov telescopes have traditionally been equipped with 
photomultiplier (PMT) detectors, a natural choice since the im¬ 
age point-spread-function in the focal plane is typically some 
cm in extent, so well matched to PMT sizes. However, their 
quantum efficiency is not perfect, and their electric requirements 
and consequences from overexposure may cause some concern. 
Recently, higher-efficiency and less fragile silicon avalanche 
photodiode arrays have been incorporated into Cherenkov cam- 


eras ( 

Anderhub et al.[[201 3[ [Catalano et al.[[20131 [Sottile et al.[ 

|20T3 

) and it can be envisioned that such detectors will also be 


used in future telescopes, at least those with an optical two- 
mirror design that produce a fairly compact focal plane. 


Single-photon-counting silicon avalanche photo-diodes 
(SPADs) have a potential for quantum efficiency that ap¬ 
proaches unity (and extends into the infrared), while counting 
individual photons at nanosecond resolution. They are oper¬ 
ated in ‘Geiger mode’ where each detected photon triggers 
an electronic avalanche as a signature of photon detection. 
Following a detection, a certain deadtime occurs, during which 
further detections are not possible. Typical deadtimes can be 
tens of nanoseconds, which permit count rates up to 10 MHz or 
somewhat higher. 

In addition to dark counts, detectors also show some level of 
afterpulsing. There remains a possibility that an avalanche elec¬ 
tron is occasionally caught in the potential well around some 
semiconductor impurity site. If that trapped electron is released 
after a time longer than the deadtime, it may trigger a new 
avalanche, which is correlated with the real photon event. Since, 
in intensity interferometry, one searches for correlated signals, 
such afterpulsing may appear as a false correlation. Of course, 
the afterpulsing is a property inside each single detector and will 
not (except as a small higher-order effect) be correlated between 
two separate ones. 

Another somewhat peculiar property of SPADs is the emis¬ 
sion of light from the detector surface (‘diode afterglow’). 
Immediately following the photon-detecting avalanche, a shower 
of photons is emitted from the detector as it recovers. For high¬ 
speed applications in large telescopes, some such secondary 
light could find its way back onto the detector, causing optical 
‘ringing’. This emission can also induce crosstalk between adja¬ 
cent pixels in a detector array (‘optical crosstalk’). In the elec¬ 
tronic design of SPAD arrays, efforts have been made to remedy 
such effects by, for instance, making trenches in the silicon to 
avoid direct exposure between adjacent pixels. 

These examples illustrate that an adequate understanding of 
detector properties is desirable in both the selection and the 
use of actual detectors. Since solid-state detectors appear to be 
a promising choice in terms of sensitivity and handling, those 
were chosen for the present experiments. One practical concern 
is that the physical size of single-pixel detectors is tiny, at only 
a fraction of a mm, so incompatible with the large image scales 
in Cherenkov telescopes, where larger array detectors have to 
be used. Large solid-state photomultiplier arrays are now avail¬ 
able from several manufacturers but already the experience from 
smaller-scale ones should be relevant for future full-scale opera¬ 
tions. 

Single-photon-counting avalanche silicon photodiode detec¬ 
tor modules (SPADs) from various manufacturers were evalu¬ 
ated in laboratory tests. Following those, a dozen modules from 
MPD, Micro Photon Devices, with sensitive areas of 100 pm di¬ 
ameter, Peltier-cooled, and with dark-count rates below 500 Hz 
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were acquired to equip a corresponding number of small labora¬ 
tory telescopes. 

3.3. Real-time photon correlation 

An essential element of an intensity interferometer is the cor¬ 
relator, which provides the temporally averaged product of the 
intensity fluctuations (A • I 2 } from two telescopes, normalized 
by the product of average intensities in each of them, (A) • (h)- 
The original interferometer at Narrabri used an analog correlator 
to multiply the photocurrents from its photomultipliers. Analog 
correl ators ( [LeBohec et al.[|2008||Nunezl|2Q12||Pellizzari et aL| 
|2012| ) have attractive properties in often being able to handle 
wide electrical band widths. However, they may also be sensitive 
to electromagnetic interference from all sorts of electrical de¬ 
vices omnipresent also in observatory environments ( [PellizzarT 
|et aH 2Q12[ ). For applications in intensity interferometry, these 
may be critical since the measured signal is expected to be tiny, 
and has to be segregated from a very much greater background. 
Of course, isolation from interference is feasible with sufficient 
shielding, and having optical cables instead of metallic ones may 
alleviate the problems, but there are still issues at their inter¬ 
faces. Because of these considerations, the signal handling was 
chosen here to only be digital. This may introduce limitations in 
how high count rates that can be handled per channel but - since 
the S/N in principle is independent of the optical bandpass (see 
below) - stellar fluxes may be adjusted to manageable photon- 
count levels using suitably narrow-band optical Alters. 

Photon correlators are commercially available for primary 
applications in light scattering from laboratory specimens 
( Becker[[2005| ). Based upon experiences from a series of hard¬ 
ware correlators from different manufacturers, current units were 
custom-made (by Correlator.com) for applications in intensity 
interferometry. These feature sampling frequencies up to 700 
MHz, are able to handle continuous photon-count rates of more 
than 100 MHz per channel without any deadtimes, and have on¬ 
line data transfer to a host computer. Their output contains the 
cross correlation function between pairs of telescopes (as well 
as autocorrelation functions for each of them), made up of about 
a thousand points. For small delays, the sampling is made with 
the shortest timesteps of just over 1 ns, increasing in a geomet¬ 
ric progression to high values to reveal the full function up to 
long delays of even seconds or more. It is believed that their 
performance is adequate for full-scale intensity interferometry 
experiments. In particular, it should be stressed that the correla¬ 
tors required here are very much more modest than those used 
in radio phase interferometer facilities (of supercomputer class), 
which decipher not only spatial but also spectral information, 
and with many bits of resolution. 

Still, such correlators are not without problems. In the ideal 
case of a spatially coherent source measured with perfect tempo¬ 
ral resolution, the average of the instantaneously squared inten¬ 
sity equals twice the square of its average (Eq.l). To obtain such 
a value in practice requires measurements with a time resolution 
shorter than the optical coherence time. While this can be real¬ 
ized for special laboratory light sources (as discussed below), it 
is not realistic for astronomical ones (except perhaps some hy¬ 
pothetical natural lasers), where the optical passband is deflned 
by some wavelength Alter. Realistic bandpasses of 1 nm or 10 
nm correspond to optical frequency bandpasses on the order of 
1 or 10 PHz and coherence times around one picosecond or less. 
Since realistic electronic resolutions are close to a nanosecond, 
the intrinsic intensity fluctuations are averaged over many opti¬ 
cal coherence times, and the measured squared intensity devi¬ 


ates only slightly from the square of its average, rather than by 
the full factor of two. With random fluctuations decreasing as 
the square root of the number of samples, this smearing may be 
a factor of 100 or more, requiring both large collecting areas to 
reach adequate photon statistics, and analysis of the correlation 
signal to many decimal digits. 

Firmware correlators produce correlation functions in real 
time, processing very large amounts of photon-count data and 
eliminating the need for their further handling and storage: for 
example, ten input channels, each running at 50 MHz during 
one eight-hour observing night, process more than 10 terabytes 
of photon-count data). A disadvantage is that, if something needs 
to be checked afterward, the original data are no longer available, 
and alternative signal processing cannot be applied. 

An alternative approach is to time-tag each photon count 
and store all data (at least for moderate count rates), and then 
later perform the correlation analyses off-line. The data streams 
from multiple telescopes can then be cross-correlated in soft¬ 
ware, possibly applying noise Alters and also computing other 
spatio-temporal parameters such as higher-order correlations be¬ 
tween three or more telescopes, which could contain additional 
information. Such a capability was built into the AquEYE and 
IquEYE instruments developed at the University of Padova for 
very high time-resolution astrophysics, also with a view toward 
inten sity interferometry, using similar SPAD detecto rs to those 
here ( [Capraro et al.[|2010}|Naletto et al.[|200^|2Q13| ). 


3.4. Image reconstruction from second-order coherence 

While intensity interferometry possesses the advantage of not 
being sensitive to phase errors in the optical light path, ordi¬ 
nary two-telescope correlations also do not permit such phases 
of the complex coherence to be measured. These correlations 
provide the absolute magnitudes of the respective Fourier trans¬ 
form components of the source image, while the phases are not 
directly obtained. Such quantitites can be used well by them¬ 
selves to fit model parameters such as stellar diameters, stellar 
limb darkening, binary separations, and circumstellar disk thick¬ 
nesses, but actual images cannot be directly computed through a 
simple inverse Fourier transform. 

While a two-component interferometer (such as the classi¬ 
cal one at Narrabri) offers only very limited coverage of the 
Fourier (u,v)-plane, a multicomponent system provides numer¬ 
ous baselines and an extensive coverage of the interferometric 
plane. Already intuitively, it is clear that the information con¬ 
tained there must place stringent constraints on the source image. 
For instance, viewing the familiar Airy diffraction pattern (cf. 
Figure 1 left), one immediately recognizes it as originating in a 
circular aperture, although only intensities are seen. However, 
it is also obvious that a reasonably complete coverage of the 
diffraction image is required to convincingly identify a circular 
aperture as the source. 

Various techniques (most unrelated to astronomy) have been 
developed for recovering the phase of a complex function when 
only its magnitude is known. Methods specifically for intensity 
interferometry were worked out by Holmes et al. (2004; 2010; 
2013) for one and two dimensions, respectively. Once a suffi¬ 
cient coverage of the Fourier plane is available, phase recov- 
ery and imaging indeed become possible. [Nunez et ^ ( 2012a|b l 
applied this phase recovery to reconstruct images from simu¬ 
lated intensity interferometry observations, demonstrating that 
also rather complex images can be reconstructed. (However, a 
limitation that still remains is the non-uniqueness between the 
image and its mirrored reflection.) 
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Fig. 1. Diffraction patterns in A 532 nm laser light show the 
[squared] Fourier transforms of some of the artihcial ‘stars’. Left 
to right: Circular single star of 150 jum nominal diameter; ellipti¬ 
cal small single star; symmetric binary star with each component 
100 jum. The side of each image corresponds to ~ 70 cm in the 
telescope plane and corresponding baselines are required to re¬ 
trieve these patterns. 


4, Telescope array in the laboratory 

The project aim was to realize an end-to-end simulation of an 
intensity interferometer observing star-like sources with a large 
array of many telescopes, equipped with electronics of the type 
foreseen for full-scale operations. An interferometer was thus set 
up in a very large optics laboratory, comprising an artihcial star 
viewed by a linear array of typically hve stationary telescopes of 
different separations. A two-dimensional telescope layout with 
many tens of telescopes could be simulated by successively ro¬ 
tating the position angle of the star relative to the plane of the 
telescopes. With high-speed photon-counting detectors and real¬ 
time digital cross correlation between many tens or even hun¬ 
dreds of telescope pairs, operations resemble foreseen future ob¬ 
servations with large Cherenkov telescope arrays. Still, such a 
laboratory setup demands several considerations. 

4.1. Small sources and long distances 

A hrst issue concerns the artihcial ‘stars’. These should have 
somewhat realistic angular sizes and also possess some non¬ 
trivial properties whether single (round or elliptical), or dou¬ 
ble (with equal or unequal components). They were prepared as 
small apertures, drilled as physical holes in metal (using minia¬ 
ture drills otherwise intended for mechanical watches). This sets 
some minimum practical size of aperture openings of ~ 100 jum. 
Such an opening observed at 20 m distance subtends an angle 
of 5-10“^ rad (1 arcsec). With the Airy disk diffraction radius 
given by 0 1.22A/D rad, the telescope aperture D required 

to resolve it comes out to ~10 cm, dictating a rather compact 
setup of telescopes unless extremely large laboratory spaces are 
available. Our interferometer occupied most of the wall-to-wall 
extent of the large optics laboratory at Lund Observatory, with 
a source-to-telescope distance of ~23 m. This choice followed 
an evaluation of alternative setups simulated with optical design 
software (ZEMAX), where we found that the requirements of 
source surface brightness and small optical aberrations to match 
a small detector size were too awkward to satisfy with variants 
involving, for example, multiple-reflection passes, optical fibers, 
or large-size lenses. 

Figure 1 shows photographs of the diffraction patterns of 
some of the artificial stars, here illuminated by coherent laser 
A 532 nm light, projected onto a large screen. These patterns 
Fourier transforms of the aperture openings serve both to ob¬ 
tain precisely the dimensions and shapes of the ‘stars’ (quantities 
otherwise difhcult to directly measure on the tiny openings) and 


to illustrate the spatial coherence patterns that can eventually be 
retrieved from interferometry. 


4.2. Dependence on source temperature 


The signal-to-noise ratio (S/N) in measuring the second-order 
coherence for one pair of telescopes equals 


(S/NhMS =A-a-n- \yn(r)f ■ (2) 


e.g., Hanbury Brown| ( [l974| ). Here A is the geometric mean of 


the areas of the two telescopes; a is the quantum efhciency of 
the optics plus detector system; rj is the flux of the source in 
photons per unit optical bandwidth, per unit area, and per unit 
time; l 7 i 2 (r)P is the second-order coherence of the source for 
the baseline vector r, with jnix) being the mutual degree of co¬ 
herence. And A/ is the electronic bandwidth of the detector plus 
signal-handling system, and T the integration time. 

Most of these parameters are easily understood since they 
directly depend on the instrumentation, but rj depends on the 
source itself, because it is a function of its brightness tempera¬ 
ture. Thus, for a given number of photons detected per unit area 
and unit time, the S/N is better for sources where those pho¬ 
tons are squeezed into a narrower optical band. This implies a 
valuable property for intensity interferometry in that the S/N is 
independent of the width of the optical passband AT, as shown 
by Eq. (2). While a smaller flux over a narrower passband wors¬ 
ens the photometric precision due to photon statistics (assuming 
a flat source spectrum), that is counterbalanced by a longer op¬ 
tical coherence time, and a diminished smearing of the physical 
signal during the electronic time resolution. 

This property has already been exploited in th e original 
Narrabri interferometer ( [Hanbury Brown et'ari|197Q| ) to identify 
the extended emission-line volume from the stellar wind around 
the WolfRayet star Vel. They compared observations in the 
A 465 nm C iii-iv emission feature emanating from the wind, with 
those in the photospheric continuum around A 443 nm, flnding 
the extent of the emission-line region several times greater than 
the stellar disk, and concluding that the stellar wind Alls the en¬ 
tire Roche equipotential lobe around this binary star. Their con¬ 
tinuum spectral region (deflned by a 10 nm wide color Alter) was 
selected to be free from prominent spectral features, while the 
ionized carbon emission was measured over a narrower band of 
2.5 nm. The same effect could also be exploited to increase the 
S/N by observing the same source simultaneously in multiple 
spectral channels, a concept foreseen for the o nce-proposed suc- 
cessor to the original Narrabri interferometer ([Hanbury Brownj 

[T^ . 


That the limiting S/N is set by the source temperature rather 
than by its brightness or telescope size, may appear somewhat 
counter-intuitive (and has been the cause of some misunder¬ 
standings) but can be interpreted in both classical and quantum 
terms. Once the observational equipment is in place, one may 
try to improve the S/N. One could try increasing the photon 
flux by going to a broader wavelength interval in white light. 
However, the S/N will not change: realistic electronic resolu¬ 
tions (~ns) are always very much slower than the temporal co¬ 
herence time of broad-band light (perhaps lO'^'^s). While broad¬ 
ening the spectral passband does increase the photon count rate, 
it also decreases the temporal coherence by the same factor. The 
intensity fluctuations have their full amplitude over one coher¬ 
ence time, and now get averaged over many more, canceling 
the effects of decreased photon noise. Then one could try with 
larger telescopes. However, once telescope sizes begin to ap¬ 
proach the structures in the spatial coherence pattern, although 
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increasing telescope areas gives more photons, they also average 
over more spatial coherence structures by the same factor, again 
not improving the S/N. Alternatively, one might be tempted to 
observe brighter sources. However, for any given source temper¬ 
ature, brighter sources will be larger in angular extent and have 
smaller structures in the Fourier domain, so one will again aver¬ 
age over additional spatial coherence areas, without improving 
the S/N. Indeed, even pointing the telescope at the Sun will not 
help since the only way to get better S/N is to find sources with 
a higher surface brightness. Below some effective source tem¬ 
perature, no sensible measurements can be achieved, no matter 
how bright the source, or how large the telescopes. A diagram 


illustrating this steep temperature dependence is in [Dravins & 
[Lagadec (2014); their Figure 3. 


4.3. High-temperature laboratory light sources 


That intensity interferometry is sensitive to sources of high 
brightness temperature but limited in observing cool ones, of 
course, is equally valid for any laboratory setup, as for stars in 
the sky. The source must be small enough to produce an ex¬ 
tended Fourier pattern that can be sampled by the interferom¬ 
eter, while also being bright enough to produce acceptable pho¬ 
ton count rates. However, while there are many stars in the sky 
with Teff = 10,000 K or more, to produce a correspondingly bril¬ 
liant laboratory source is much more challenging. Furthermore, 
the precise relation between second- and first-order coherence, 
= 1 + assumes that the light is chaotic, with a 

Gaussian amplitude distribution ( |Shih[ [201 1[ [Bachor & Ralphj 
|2004t Loudon|[2000| ), i.e., the light waves undergo random phase 
shifts so that intensity fluctuations result, which then bear a sim¬ 
ple relation to the ordinary first-order coherence. While this must 
be closely satisfied for any thermal source, it is not the case for 
lasers that, ideally, never undergo any intensity fluctuations any¬ 
where = 1). The spatial extent of a laser source therefore 
cannot be measured by intensity interferometry, and a laser is 
not an option for enhancing the brightness of an artificial star. 


In initial attempts to achieve a high surface brightness for 
an artificial star, the very small emission volume of a high- 
pressure Hg arc lamp was focused onto a pinhole serving as 
the star, and a narrow-band optical filter singled out its bright¬ 
est emission line (A 546 nm). Although this represents about 
the highest blackbody brightness temperature (some 5,500 K) 
that can be obtained with somewhat ordinary laboratory equip¬ 
ment for a non-laser source, the photon-count rates still turned 
out to be too low for conveniently short integration times. Other 
tests were made with several lamps of various atomic species, 
selecting their brightest and narrowest emission lines, but the 
signals still remained marginal. Finally, quasi-monochromatic 
chaotic light was produced by scattering monochromatic laser 
light from microscopic particles, suspended in a water-filled cu¬ 
vette. Such particles undergo random thermal (Brownian) mo¬ 
tion due to collisions with the water molecules that surround 
them, producing a slightly Doppler-broadened spectral line that 
is extremely narrow (~ MHz). Thus, a light source with an ef¬ 
fective brightness temperature estimated as Teff -100 000 K was 
produced, permitting conveniently short integration times on the 
order of a minute. Such scattered laser light is used for various 
laboratory photon-correlation measurements of time variability 
( [Becker I [2005[ ), but we are not aware of any previous such ex¬ 
periment with a spatial intensity interferometer. 



Fig. 2. Measured second-order temporal coherence of light used 
to illuminate the artificial ‘stars’. Laser light was scattered 
from suspensions of plastic microspheres undergoing thermal 
(Brownian) motion in water. The coherence falls to unity for 
delays longer than the characteristic coherence times. These 
are longer for larger particles (curves with longer dashes), 
which move more slowly, inducing less Doppler broadening. 
Analogous functions are seen for diluted milk (solid curve with 
all measured data points), where microscopic fat globules act 
as scatterers. These measurements of temporal coherence were 
made with a single telescope on a small and spatially coherent 
source. 


4.4. Dynamic light scattering 


Monochromatic light scattered by microscopic particles in ther¬ 
mal motion is broadened by Doppler shifts, producing a spectral 
line with a Lorentzian wavelength shape, and with a Gaussian 
(thermal, chaotic, maximum-entropy) distribution of the elec- 
tric fi eld amplitudes ( [Berne & Pecora[ [2000} [Crosignani et al. 


1975[ ). This light is thus equivalent in its photon statistics anc 
intensity fluctuations to the thermal (white) light expected from 
any normal star, the difference being that the spectral passband 
is now very much narrower. In addition to a high Teff, the nar¬ 
rowness of the spectral line implies a long coherence time, less¬ 
ening the requirements on the detector time resolution and per¬ 
mitting measurements also on (the less demanding) microsec¬ 
ond scales, where detector concerns such as afterpulsing are less 
pronounced. The exact values for the spectral broadening can be 
modified by choosing the particle size (larger particles undergo 
slower motion, inducing less Doppler broadening), but it also 
somewhat depends on the viscosity and the temperature of the 
medium of suspension. 

Such light sources are not completely without problems, 
however. Intensity interferometry relies upon a relation be¬ 
tween first- and second-order coherence that is valid for chaotic 
light. The scattering from single and non-interacting particles 
in Brownian motion produces this light, but only in the case of 
single scattering. In the cases where the photons undergo mul¬ 
tiple scattering events, if the scattering particles experience mu¬ 
tual interactions, where they have intrinsic dynamic properties, 
or where there are internal currents in the suspension liquid, this 
relation may no longer precisely hold. In such cases, there will 
still be relations between and g^^\ but possibly in a some¬ 
what different form. To avoid such complications, the optical 
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Fig. 3. Components in the laboratory setup. Left: Light from a 300 mW A 532 nm laser is made chaotic by scattering from mi¬ 
croscopic particles in a square-top cuvette and focused by a condenser onto artificial ‘stars’, which are mechanical apertures in a 
rotatable holder. Right: The ‘stars’ are observed by a group of (here) five small telescopes with 25 mm apertures, each equipped 
with a photon-counting SPAD detector. One unit perpendicular to the others uses a 45-degree mirror to obtain a particularly short 
baseline. Another pair of two telescopes behind one beamsplitter serves to measure zero baseline, as required for calibrations. Two- 
dimensional coverage is achieved by successively rotating the position angle of the source relative to the plane of the telescopes. 


paths were minimized by focusing the laser light very close to 
the surface of the fiuid inside the cuvette. 

Chaotic light was thus produced by scattering T 532 nm laser 
light from microscopic monodispersive polystyrene spheres 
(from Polysciences) suspended in a cm-sized cuvette with room- 
temperature distilled water. After trying out spheres of several 
different sizes (ranging between 50 nm and 90 jim), a diameter 
of 0.2 jum was chosen as giving both an acceptably high light 
level in the scattering, and a convenient coherence time for mea¬ 
surements (Figure 2). 

However, the experiments required quite a number of tests, 
where the cost of volumes of such microscopic spheres started to 
become a problem. As an inexpensive scattering fiuid with com¬ 
parable properties, homogenized household milk (diluted with 
two parts water) was then used instead. Milk contains micro¬ 
scopic fat globules, which undergo analogous Brownian motion. 


also studied with photon correlation spectroscopy ( 

Alexander &| 

Dalgleish| 2006[ Holt et al.||1973[ Robin & Paquin| 

\99\). Given 


that the cost per volume is orders of magnitude lower than that 
of microsphere suspensions, several of the initial experiments 
were carried out with laser light scattering in milk. Also, milk 
with varying fat content was tried out, but with no obvious dif¬ 
ferences in terms of, say, correlation timescales. Presumably, the 
homogenization process produces fat globules of similar size, 
and the difference between milk with differing fat content lies in 
its quantity, not its microscopic properties. 

The final optical setup was preceded by several experi¬ 
ments and optical ray-tracing simulations. Only a relatively nar¬ 
row range of parameters was found, where laboratory measure¬ 
ments appeared practical, in particular constrained by achievable 
source brilliance. The laser was focused by optics into a very 
small volume in the scattering fiuid, and that small illuminated 
volume was focused by a large condenser onto the tiny aperture 
of the artificial star, its brightness optimized with a three-axe 
micrometer. For example, smaller artificial stars would have per¬ 
mitted a less compact and more easily handled interferometer ar¬ 
ray. However, decreasing the pinhole diameter by just a factor 2, 
say, not only decreases the area and photon flux for the artificial 
star by a factor 4, but also diffracts its light over a twice greater 


angle, thus four times greater area, leading to an illumination 
decrease in the telescope plane by a factor 16. Since the signal 
for intensity correlation is proportional to intensity squared, the 
observed signal would decrease by a factor 256, leading to im¬ 
practical integration times. 

The chosen telescopes were a compromise between larger 
sizes to collect more photons and smaller ones to fully resolve 
the Fourier-plane structures. Again, the functional dependences 
are steep: decreasing the telescope diameter a factor of 2, de¬ 
creases its light-collecting area a factor 4 and the correlation sig¬ 
nal by 16. With careful optical design and using a more power¬ 
ful laser, one should be able to increase the source brilliance, but 
probably not by a very large factor: focusing many watts into 
the submillimeter volume of the scattering fiuid might induce 
currents, clog the microspheres, induce bubbles or even cause 
boiling. 


4.5. Interferometer setup 


Figure 3 shows some photographs of the laboratory setup. The 
telescopes are small refractors with 25 mm diameter achromatic 
objective lenses mounted on an optical bench at some 23 me¬ 
ters distance from the artificial stars. In each telescope, light is 
focused onto a SPAD, a single-photon-counting avalanche pho¬ 
todiode, enabling photon-count rates up to ~10 MHz. The pulse- 
train output (electronic TTL standard) is fed to a computer- 
controlled firmware correlator for real-time cross computation of 
the data st reams between various pair s of telescopes. For more 
details, see [Dravins & Lagadec|(|2014|). 


4.6. Measurements of single ‘stars’ 

Examples of measured intensity cross correlation functions be¬ 
tween a pair of telescopes, (hit) • hit At)}, are in Figure 4, 
normalized to = 2 for zero telescope baseline. For long delay 
times At, the correlation vanishes, and the values tend to unity, 
as is appropriate for random uncorrelated variations. The char¬ 
acteristic slope for delays around Ar ~10 yus is a measure of the 
temporal coherence time and indicates the optical bandpass of 
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Fig. 4. Cross correlation functions measured for an artificial sin¬ 
gle star of 1 arcsec apparent diameter (100 /im aperture at 23 m), 
normalized to zero baseline, show how (a) the temporal coher¬ 
ence gradually decreases with increasing time delay, and (b) the 
spatial coherence gradually decreases with increasing baseline. 
In this particular case, the normalized value for zero baseline was 
taken as the autocorrelation for delays in the range 1-10 yus (at 
much shorter delay times the autocorrelation rises to unphysical 
values due to detector afterpulsing). These effects were avoided 
in other measurements where the zero baseline was realized by 
cross correlating two telescopes behind one beamsplitter. 



Interferometer baseline [cm] 


Fig. 5. Second-order spatial coherence measured for artifi¬ 
cial single stars of different angular sizes. The coherence val¬ 
ues come from correlation functions such as in Figure 4, aver¬ 
aged for delays between 1-10 yus, normalized to baseline zero. 
Superposed are theoretical dependences for ideal circular aper¬ 
tures (Airy patterns; squared moduli of the Fourier transforms, 
1+ (2 Ji(x)/x)^), confirming the expected trends. The precise 
effectively illuminated source shapes depend on focusing details 
and are not exactly known; however typical measuring repro¬ 
ducibilities were within ~0.02 of the coherence values. 


the scattered laser light, here on the order of 10^ Hz. For broader 
band light, the coherence time is much shorter, and then the pos¬ 
sibilities of resolving the varying correlation for successively 
shorter delays become limited. However, in this arrangement it is 
possible to follow the detailed changes in both spatial and tem¬ 
poral correlations, when the temporal averaging proceeds over 
successively larger numbers of coherence times. 


For intensity interferometry we are not concerned about the 
temporal coherence but only the spatial one: its gradual change 
with increasing baseline enables the angular size of the source 
to be determined. At delays shorter than ~1 yus, the temporal 
coherence is here fully resolved, and the differences between 
curves then reflect the spatial coherence only. For zero base¬ 
line, the spatial coherence must approach = 2, and the auto¬ 
correlation function is normalized accordingly since several ef¬ 
fects contribute to the actually measured values being lower. For 
example, a slight deviation from zero baseline occurs because 
the telescope dimensions are not completely negligible and their 
measured spatial coherence corresponds to the convol ution of 
the telesc ope pupils with the spatial coherence pattern ( Rou et| 


al. 2013| . The ideal value of g^^^ = 2 is retrieved only for fully 


polarized light, while in the u npolarized light used he re, the sig¬ 
nal decreases by a factor ^/2 (Hanbury Brown 1974| . Also, the 
measured signal is diluted by some (small) fraction of coherent 
laser light that did not scatter from the suspended microparticles 
but illuminated the aperture by being refracted in the walls of 
the glass cuvette, preserving its second-order coherence g^^^ = 
1. To account for these effects, the successive cross correlation 
functions were normalized relative to the autocorrelation ones. 
Figure 5 shows the dependences of the measured second-order 
coherence on telescopic baseline. 


4.7. Hundreds of baselines: Mimicking large arrays 


A number of more complex sources measured included asym¬ 
metric stars and binary sources of varying sizes and separa¬ 
tions. To mimic large telescope arrays, the effective number of 
telescopes and baselines was increased by making successive 
series of measurements across different angles in the interfer¬ 
ometric (u,v)-plane. The telescope array is mounted horizon¬ 
tally and thus only covers horizontal baselines. To effectively 
cover oblique and vertical baselines as well, the artificial star 
was rotated to successively different position angles to change 
the (u,v)-plane coverage. With N telescopes at non-redundant 
mutual separations, N(N-l)/2 different baselines can be con¬ 
structed. With five telescopes, ten baselines are available at each 
angular position and, for example, ten different position angles 
produce 100 different baselines. The second-order coherence in 
a point (u,v) of the Fourier-transform plane equals the one in 
(-u,-v), so each baseline provides two points. Figure 6 shows 
coherence patterns for sources measured over 60 and 100 base¬ 
lines, and Figure 7 such from 180 optical baselines. 

With such numbers, one moves into new parameter domains 
for optical interferometry, beginning to approach the capabili¬ 
ties of large Cherenkov arrays. Assuming 50 telescopes will be 
available for interferometry, the number of baselines, N(N-l)/2, 
would already exceed 1000 (even if some might be redundant 
due to repetitive telescope locations). As verified in numerical 
simulations, the ensuing nearly complete (u,v)-plane coverage 
enables full two-dimensional image reconstruction (|Nunez et 
aL] |2012a|bp . One additional requirement during astronomical 
observing (although not needed here) will be to electronically 
and continuously track the changing projected baselines between 
pairs of telescopes as the source moves overhead across the sky 
during an observing night, assigning measurements to the instan¬ 
taneous location in the (u,v)-plane. 
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sity interferometry has the advantage of not being sensitive to 
phase errors in the optical light paths, it also does not measure 
such phases but instead obtains the absolute magnitudes of the 
respective Fourier-transform components of the source image. 
These data can fit model parameters such as stellar diameters 
or binary separations but images cannot be computed through 
merely an inverse Fourier transform. 

Some strategies that specialize in anal yzing intensity inter- 
ferometry data have been developed; e.g., Holmes & Belen’kii 
( |2004[ ); Nunez et ‘^^ ( [2012^ , and consist in estimating Fourier 
phases from Fourier magnitudes. Since the Fourier transform of 
a finite object is an analytic function, the Cauchy-Riemann equa¬ 
tions can be used to find derivatives of the phase from deriva¬ 
tives of the magnitude, producing images that are unique ex¬ 
cept for translation and reflection. Here, one could start with a 
one-dimensional phase estimate along a single slice through the 
(u,v)-plane origin. Two-dimensional coverage requires combin¬ 
ing multiple 1-D reconstructions, while ensuring mutual consis¬ 
tency between adjacent ones. However, these algorithms require 
a very dense coverage of the (u,v)-plane and - at least on first 
attempts - did not prove efficient with current laboratory data. 

Therefore, we chose to follow an inverse-problem approach, 
where a priori constraints are imposed to interpolate between 
successive Fourier frequencies. The imposed constraints are very 
general in character, and do not require much knowledge of the 
object. Such prior constraints (known as ‘regularization’) are 
implemented i n the ‘multi-aper ture image reconstruction algo¬ 
rithm’, MiRA: |Thiebaut| ( |2009| ). It has be en tested for analyz - 
ing simulated intensity interferometry data ( [Nunez etaT] 2012b| ), 
and we use it here. The MiRA maximizes agreement with data, 
as well as with a general penalty (regularization) function. It can 
be stated as the following constrained optimization problem for 
an image x: 


'X-best = argmirijc {fdatuix) + flfpriorix.)] 


(3) 


Fourier u = k^X [cm] 

Fig. 6. Second-order spatial-coherence surfaces measured at 
T 532 nm for artificial stars in the laboratory. Top: Binary ‘star’ 
with equal components of diameter ~1 arcsec, from intensity 
correlations measured between pairs of telescopes across 60 dif¬ 
ferent non-redundant baselines. The pattern of central maxima 
indicates the binary separation while the symmetric rings re¬ 
veal the size of individual stars. Black dots mark individual data 
points, also used in the later image reconstructions (the surfaces 
merely connect these measured points). Bottom: A somewhat ir¬ 
regular and elliptical ‘star’ of just below 1 arcsecond, measured 
over 100 baselines. The data span the interferometric (u,v)-plane 
of the Fourier transform of the source image with different spa¬ 
tial (angular) wavenumbers k, corresponding to different tele¬ 
scope baseline lengths and orientations. These maps provide the 
[moduli of the] Fourier transforms of the source images, analo¬ 
gous to the first-order diffraction patterns of intensity produced 
in coherent light (Figure 1). Images reconstructed from these 
data are in Figure 8. 


5. Image reconstruction 

Image reconstruction from second-order coherence implies 
some challenges not present in first-order phase/amplitude tech¬ 
niques, and remains an active subject of research. While inten¬ 


such that > 0 for all k and 

Here x is the image, Xk the pixel value, and yu a constant. 
The function f^ata quantifies the agreement with the data as 
a while /prior is the regularization function. Here we use 
two different regularizations: smoothness and compactness, in¬ 
cluded in MiRA. The smoothness penalty function is defined 
as /prior = |x - S.xp, where S is a smoothing operator imple¬ 
mented via finite differences. The compactness penalty is de¬ 
fined as /prior = The variable quantifies the global weight 
of the penalty function. The additional constraints are the pos¬ 
itivity and normalizability of the image. For reconstructing the 
image of the single elliptical star from the coherence patterns in 
Figure 6 (bottom), the compactness regularization was used with 
yu = for the symmetric binary (Figure 6, top) the smooth¬ 
ness regularization was applied with yu = 10^^, and for the asym¬ 
metric one (Figure 7) yu = 10^^. The resulting reconstructed im- 
ages are in Figure 8, for which a shorter description appeared in 
[Dravins et al.| ( |2015] ). 

Besides limits set by measurement noise, the fidelity of any 
reconstruction depends on how completely the (u,v)-plane has 
been filled with data. In practice, every spatial frequency cannot 
be measured, and a fit to any sparse data set cannot yield a com¬ 
pletely unique image without artifacts. Of course, this also holds 
for any image reconstruction from a finite number of baselines 
in ordinary amplitude/phase interferometry. 


9 

















































































































D. Dravins et al.: Long-baseline optical intensity interferometry 



Fig. 7. Second-order optical coherence of an artificial asymmet¬ 
ric binary star with differently sized components, built up from 
intensity-correlation measurements over 180 baselines between 
pairs of small optical telescopes. The coordinates refer to the 
plane of the Fourier transform of the source image and corre¬ 
spond to different telescope baseline lengths and orientations. 
At bottom, the projection of the 3-D mesh is oriented straight 
down, showing [the modulus of] the sources Fourier transform 
(diffraction pattern). The image reconstructed from these data is 
in Figure 8, bottom. 


5.1. Image reconstruction practicalities 

The numerical stability and reproducibility of the reconstruc¬ 
tions were examined. For the single asymmetric star in Figure 
8 (top), its semi-major axis was measured as 0.53 + 0.01, and 
its semi-minor one as 0.36 + 0.01 arcsec. The ‘edge’ of the star 
was here taken as pixels above 0.5% of the image maximum, and 
the reproducibility estimated by initiating the image reconstruc¬ 
tion with different random images. However, the current setup 


10 



Fig. 8. Images reconstructed from intensity interferometry mea¬ 
surements of (top to bottom) a single, somewhat irregular el¬ 
liptical ‘star’, a binary ‘star’ with equal, and one with unequal 
components, from measurements with 100, 60, and 180 optical 
baselines, respectively. Sufficiently dense coverage of the inter¬ 
ferometric Fourier-transform plane enables a reconstruction of 
the full source images, despite the lack of Fourier phase infor¬ 
mation. The scale on the right denotes relative intensities. In 
principle, the images are unique except for their possible mir¬ 
rored reflections. These are believed to be the first diffraction- 
limited images that have been reconstructed from an array of 
optical telescopes, connected only through electronic software, 
with no optical links between them. 
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does not permit a very precise comparison of the retrieved im¬ 
age to the original since the effective source shape is set by the 
illumination pattern falling onto the geometric aperture. While 
the latter was independently measured in homogeneous laser il¬ 
lumination (Figure 1), the critical focusing of its illuminating 
scattered light probably produces some inhomogeneities across 
it. 

In the case of the binary stars, one limitation comes from the 
finite (u,v)-plane coverage: with baselines between 3 and 20.5 
cm, the former constrains the largest field to nominally 4.5 arc- 
sec, and the latter limits the resolution to 0.65 arcsec. Precise 
measures that can still be obtained include the angular separation 
of the symmetric binary, measured between the highest-valued 
pixels as 1.90 + 0.01 arcsec. For the asymmetric binary, its flux 
ratio 4) is retrieved reasonbly well, but the relative sizes of 
the components are not. This behavior is consistent with earlier 
simulations, which found that the precision is affected by the 
brightness ratio between both binary members, becoming an is¬ 
sue when one of them is more than about three times brighter 
than the other ( [Nunez et'SH|2012a'| ). 

The computing effort is modest. Starting from a random 
image, and measurements across 100 baselines, a 10 , 000 -pixel 
image is retrieved in some 10 seconds on a 3 GHz processor. 
This time can be shorter if there already is some initial guess 
that somewhat resembles the final image (perhaps obtained from 
some ‘raw reconstruction from a Cauchy-Riemann analysis), or 
can be longer for additional baselines or more numerous pixels. 
In any case, it is by far not comparable to the extensive numeri¬ 
cal efforts involved in radio aperture synthesis that also involve 
spectral analysis for each spatial point. 

The quality of the reconstructed images may depend on the 
density of numerical sampling across the Fourier plane. For the 
reconstructions in Figure 8 , no smoothing or interpolation of the 
measured data points were carried out. Tests with various sam¬ 
pling densities indicated that for best image reconstruction, the 
number of (u,v)-plane samples should roughly equal the number 
of baseline measurements. Possibly, other dependences could 
hold for noisier data, where some smoothing could be called 
for. In general, however, the optimization of reconstruction in 
the presence of various noise sources and incomplete Fourier- 
plane coverage is a fairly complex problem that requires further 
study (analogous to what in the past has been done for aperture 
synthesis in radio). 

Similar to any aperture synthesis from phase/amplitude in¬ 
terferometry, the fidelity of the reconstructed images depends on 
the density of measurements and on their extent across the (u,v)- 
plane, as well as the performance of any ‘cleaning algorithms 
applied to remove signatures of the (u,v)-plane sampling. While 
such algorithms of course cannot create information content on 
spatial scales that are not measured, properly deconvoluted im¬ 
ages will not display confusing instrumental signatures. 

5.2. Triple correlations and ‘closure Intensities’ 

Besides the cross correlation required to retrieve the spatial co¬ 
herence there is a potential to also explore higher-order cor¬ 
relations in light. Current measures are obtained as cross cor¬ 
relations between intensity fluctuations at two spatial locations, 
at one instant in time. However, using telescope arrays, one can 
construct, for example, third-order intensity correlations, for 
systems of three telescopes: (I(ri,ti)-I(r 2 , t 2 )-I(r 3 , ^ 3 )), where the 
temporal coordinates t do not even have to be equal. In principle, 
such higher-order correlations in light carry additional informa¬ 
tion about the source ( |Fontana||l9^||Gam^|1963t|Malvimat etj 


al. 2014[|Marathay et aLl|1994t[Qfir & Rib^|2006[|Sato et al." 


1978[|1981| ). For instance, from correlations among all possible 
triplets and quadruplets of telescopes, a more robust full recon¬ 
struction of the source image should be possible. Here, we are 
entering largely unexplored territory for observational astron¬ 
omy, but one that could well be enabled by the availability of 
extended arrays with large and numerous telescopes. 

A two-point correlation measures the squared modulus of 
the complex degree of mutual coherence 7ij7*j between two 
telescope locations i and j. Analogously, it can be shown that 
three-point intensity correlations allow a quantity to be mea¬ 
sured that is proportional to the product of the complex de¬ 
gree of mutu al coherence between detectors in a closed loop, 
i.e., 7ij7jk7ki (jMalvimat et al.||2014|). By expressing jij as 7 ij = 
\ 7 ij\exp{i((pij -t-Ajy)|, where 0// is the astrophysical Fourier phase, 
and Aij is the atmospherically induced phase, it is straightfor¬ 
ward to show that the atmospheric phases cancel out in the 
triple product. The phase of the triple product is a linear com¬ 
bination of the Fourier phases, and is known as the ‘closure- 
phase’ in amplitude interferometry. Measuring the closure phase 
for many triplets in an interferometric array allows the Fourier 
phase across the array to be constrained, which may in t urn en- 
hance image reconstruction. jMalvimat et aL ( 2014| ) and Wentz 
& Sahaj ( |2015| ) investigated the S/N properties of three-point 
intensity correlations, showing that this quantity may be mea¬ 
surable for bright stars, while Nunez & Domiciano de Souza 
(in preparation) are investigating the use of three-point corre¬ 
lations for enhancing image reconstruction. While three-point 
correlations may turn out to provide useful information about 
the Fourier transform of the source, these have more severe S/N 
limitations than do two-point ones. Therefore, their main role 
in imaging should be to supplement two-point data. However, 
if such closure-phase intensity interferometry turns out to be 
useful, it could provide additional constraints at very little ad¬ 
ditional effort since the intensity signals from each telescope can 
be freely copied or even stored off-line for later analysis. 


6. Toward kilometer-scale optical interferometry 


The current laboratory setup of a multitelescope array has 
demonstrated the end-to-end operation of an imaging intensity 
interferometer of the Hanbury Brown-Twiss type, involving very 
many baselines and densely populating large areas of the inter¬ 
ferometric Fourier (u,v)-plane. The photon count rates achieved 
and handled by the correlators in real time (^ 1 MHz) are compa¬ 
rable to what can be expected in actual stellar observations, with 
some spectral feature selected through a narrow bandpass filter. 
Already our current correlators can handle signals from up to 20 
telescopes simultaneously and, together with experiments car¬ 
ried out with actual Cherenk ov telescopes ( [Dravins & LeBohe'c 


2008[[LeBohec et ak) 201Q[ ), these experiences bode well for fu 


ture full-scale observations of celestial sources. 


6.1. Pending laboratory issues 

The current experiments were carried out with thermal light, 
which is produced by laser scattering. Although the photon 
statistics are believed to closely match those of thermal white 
light, the coherence times are much longer than the subnanosec¬ 
ond scales relevant for actual astronomical sources. Since co¬ 
herence times for broader-band light are very much shorter 
than realistic electronic timescales, most of the cross correla¬ 
tion signal will then be confined to the shortest temporal de¬ 
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lays of timescales on a few nanoseconds. Here, various detec¬ 
tor issues begin to appear, such as deadtimes and afterpulsing, 
who se understanding, calibratio n and mitigation may be essen¬ 
tial ( Dravins & Lagadec 2014 ). With the aim of realizing labo¬ 
ratory experiments also on such short timescales, some tests with 
high-temperature white-light lamps have already been made, in¬ 
dicating that it might be feasible to achieve meaningful S/Ns, 
although some effort will be required. 

For the observability of fainter sources, the S/N can be 
improved by simultaneous measurements in multiple spectral 
bands, exploiting the property that S/N is independent of opti¬ 
cal bandpass. Possibly, some scheme for dispersive optics could 
be devised, probably near t he focal plane , follow ed by a multi¬ 
element detector array; e.g., Gerwe et al. ( |2013| ). A better solu¬ 
tion, however, would be to obtain detectors that are both energy¬ 
resolving and photon-counting. While such devices have been 
demonstrated in the laboratory, they do not yet seem to be prac¬ 
tical for field work with large telescopes. Obviously, the highest 
quantum efficiency is desirable, such as values ^90% that have 
been demonstrated with superconducting nanowire devices for 
single-photon detection, reaching into the infrared. Also polar¬ 
ization properties could be optimized: Since the full signal 
appears in polarized light, that could be exploited by introducing 
a polarizing beamsplitter and detecting each polarization sepa¬ 
rately. 

How many of these features that can be practically com¬ 
bined will determine what limiting magnitudes ultimately can 
be reached, with a potential for going substantially fainter than 
the my=8 estimate for single-wavelength observations discussed 
in Section 3 above. Further information is contained in various 
higher-order correlations, but these have not yet been experimen¬ 
tally studied in the context of astronomical imaging, and theoret¬ 
ical simulations are lacking here also. 


6.2. Outlook for field experiments 

Cherenkov telescopes normally produce cm-size stellar images 
in their focal planes, which are much larger than the single-pixel 
SPAD detectors used here. Although our detectors are similar 
in principle to those large-area solid-sta te photomultipliers now 
demonstrated in Cherenkov telescopes dAnderhub et al.| 2013[ 


[Catalano et al.||2Q13}|Sottile et al. |2Q13| ), those larger area ones 

are built up by a large number of individual light-sensitive ar¬ 
eas on the same silicon chip and (as we have examined in other 
laboratory experiments) possess somewhat different deadtime, 
dark-count, and afterpulsing characteristics. An adequate under¬ 
standing of these detector properties (or of photomultipliers, in 
case such would be used) will be required to reach the photomet¬ 
ric precisions required for disentangling physical values of 
from stochastic or systematic noise. 

When observing celestial sources during an observing night 
with stationary telescopes, projected baselines between pairs of 
telescopes will gradually change. Of course, this enables an even 
richer (u,v)-plane coverage, and in deed it is the principle of aper¬ 
ture synthesis using Earth rotation ( |Saha|[MTT] ). However, to re¬ 
tain a high Fourier-plane resolution in real-time correlation, this 
requires an electronic unit that implements a variable time delay 
onto the stream of photon pulses, compensating for the relative 
timings of the wavefro nt at the different teles copes, as the source 
moves across the sky ( [Dravins et al.j [2012[ ). In case the corre¬ 
lation is computed off-line, following a recording of the signal 
during observation, the variable delay would instead have to be 
applied in software. 



Fig. 9. Vision of microarcsecond optical imaging: Expected res¬ 
olution for an assumed transit of a hypothetical exoplanet across 
the disk of Sirius. Stellar diameter =1.7 solar; distance = 2.6 pc; 
angular diameter = 6 mas; assumed planet of Jupiter size and 
oblateness; equatorial diameter = 350 //as; Saturn-type rings; 
four Earth-size moons. With a Cherenkov telescope array span¬ 
ning 2 km, already a 50 juas resolution provides more than 100 
pixels across the stellar diameter. 


7. Microarcsecond astrophysics 

7.1. Observing with small Cherenkov telescope arrays 

Probably the ‘easiest’ targets for intensity interferometry ob¬ 
servations are relatively bright and hot, single or binary stars 
of spectral types O and B or Wolf-Rayet stars with their vari¬ 
ous circumstellar emission-line structures. Their stellar disk di¬ 
ameters are typically -0.20.5 mas, so they lie (somewhat) be¬ 
yond what can be resolved with existing instruments. For ampli¬ 
tude interferometry, limitations also arise because of the paucity 
of unresolved (but sufficiently bright) reference stars for cali¬ 
brating fringe visibilities. Intensity interferometry does not re¬ 
quire such reference sources but merely instrumental stability 
between measurements made with shorter and longer baselines, 
very likely becoming critical in enabling baselines on the order 
of 1 km or more. 

Other candidate targets include rapidly rotating stars, with 
oblate shapes deformed by rotation, circumstellar disks, winds 
from hot stars, blue supergiants and extreme objects, such 
as T] Carinae, interacting binaries, the hotter parts of [su¬ 
per] nova_£xplosions,pi^atingCepheids^_or other vari¬ 

ables ( [Dravins et al.[[2012|[2Q13t[Trippe et al]l[2Q14[ ). 

Given that the S/N is independent of optical bandwidth, one 
might just as well observe in the light of some specific spectral 
feature, as in white light. Assuming suitable wavelength filters, 
one could try to map the non-radial pulsations across the sur¬ 
faces of stars such as Cepheids. The amplitudes in temperature 
and in white light probably are modest but the associated ve¬ 
locity fluctuations might be observable. If the telescope optics 
can be collimated enough to permit the use of very narrow-band 
spectral filters centered on strong absorption lines, the local stel¬ 
lar surface will appear at their particular residual intensity when 
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at rest relative to the observer, but will reach continuum intensity 
when local velocities have Doppler-shifted the absorption line 
outside the filter passband. Once such spatially resolved obser¬ 
vations of stellar non-radial oscillations have been realized, they 
should provide significant input to models of stellar atmospheres 
and interiors. 


7.2. Targets for large Cherenkov telescope arrays 


The full Cherenkov Telescope Array is foreseen to have its nu¬ 
merous telescopes distributed over a few square kilometers, with 


an edge-to-edge distance of two or three km (Acharya et al. 
|2013t CTA[ |2015| ). If fully equipped for intensity interferome¬ 
try at the shortest optical wavelengths, the spatial resolution will 
approach ~30 //as. Such resolutions have hitherto been reached 
only in radio, and it is problematical to speculate on what fea¬ 
tures could appear in the optical. However, to appreciate the 
meaning of these resolutions. Figure 9 shows an ‘understand¬ 
able’ type of object: a hypothetical exoplanet in transit across 
the star Sirius (Tgff = 9,940 K). The planet size and oblateness 
was taken as equal to that of Jupiter, but fitted with a Saturn-like 
ring and four Galilean moons, here assumed to be Earth-sized. 
Such an exoplanet would subtend a readily resolvable angle of 
some 350 //as. While spatially resolving the disk of an exoplanet 
in its refiected light may remain unrealistic for the time being, 
the imaging of its dark silhouette on a stellar disk while cer- 
tainly still very challenging might not be impossible (|Strekalov| 
|etaLl|20T3|[2tol ). 


Intensity interferometry actually possesses some advantages 
for these possible observations. The lack of sensitivity to the 
phases of the Fourier components of the image could be an ad¬ 
vantage here since one would always measure ‘only’ the ampli¬ 
tude of the Fourier transform of the exoplanet image, irrespec¬ 
tive of where on the stellar disk it happens to be. Unless close 
to the stellar limb, the star only serves as a bright background, 
its comparatively ‘huge’ diameter of 6 mas not contributing any 
sensible spatial power at any relevant telescopic baselines. 
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